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I We investigate the effects of colour-octet contributions to the radiative T decay within 

the Bodwin, Braaten and Lepage NRQCD factorization framework. Photons coming both 
^H. from the coupling to hard processes ('direct') and by collinear emission from light quarks 
r~| ! ('fragmentation') are consistently included at next-to-leading order (NLO) in a^. An estimate 
' for the non-perturbative matrix elements which enter in the final result is then obtained. By 
■ '-j ■ comparing the NRQCD prediction at NLO for total decay rates with the experimental data, 
[ it is found that the non-perturbative parameters must be smaller than expected from the 
' naive scaling rules of NRQCD. Nevertheless, colour-octet contributions to the shape of the 
photon spectrum turn out to be significant. 
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1 Introduction 



Since the early times of QCD, heavy quarkonia decays have been considered among the most 
promising processes to test the perturbative sector of the theory and to extract the value of the 
strong coupling at scales of the order of the heavy-quark mass. In addition to the calculation 
and comparison of full inclusive decay rates, much attention has been devoted to the decays 
in which one photon is emitted, and its energy measured Experimental data on the direct 
photon spectrum in T decays have been compared p|, |^, up to now, under the assumption of a 
factorization between a short-distance part describing the annihilation of the heavy-quark pair 
in a colour- singlet state and a non-perturbative long-distance factor, related to the value of the 
non-relativistic wave function at the origin. 

Recently Bodwin, Braaten and Lepage (BBL) [^] provided a new framework to study quarko- 
nium decay within QCD. Introducing an effective non-relativistic theory (NRQCD), perturbative 
and gauge-invariant factorization is obtained by including in the decay intermediate QQ states 
with quantum numbers different from those of the physical quarkonium state. The relative im- 
portance of various contributions depends on short-distance coefficients which are calculable by 
standard perturbative techniques, and on long-distance matrix elements, which can be either ex- 
tracted phenomenologically from the data or calculated on the lattice. In the end one is able to 
organize all these terms in a double perturbative series in the strong coupling and in the relative 
velocity v of the heavy quarks, and then to make predictions at any given order of accuracy. 

In quarkonia decays, photons arise from electromagnetic coupling to both heavy and light 
quarks. While contributions coming from the former, at leading order (LO) in as in the Colour- 
Singlet Model (CSM) -i.e. at the lowest order in v expansion in NQRCD- have been known 
for a long time and are one of the first tests of QCD H, LO contributions coming from 
coUinear emission from light quarks have surprisingly been considered only recently by Catani 
and Hautmann [|5|. The inclusion of these 'fragmentation' contributions within the CSM was 
found to greatly affect the photon spectrum in the T decay at low values of the energy fraction 
taken away by the photon |^. Moreover, one finds that at LO such a contribution comes entirely 
from the gluon, as the decay into light quarks vanishes. 

It then becomes natural to assess to which extent this picture remains unchanged at next- 
to-leading order (NLO) in as and in v. The aim of this work is to investigate the effects of 
colour-octet intermediate states on the photon spectrum, at fixed order in a^, aem, "f^, including 
the coupling of the photons to light quarks and gluons. In fact, while the order of magnitude of 
octet contributions is predicted using scaling rules, and found to be suppressed by powers of v 
with respect to the LO colour-singlet ones, their short- distance coefficients receive contributions 
at lower order of as-, and are then numerically enhanced. Furthermore, once leading logarithmic 
corrections are included, it is found that, contrary to the colour-singlet case, quark and gluon 
fragmentation into a photon appears at the same order in the a^, ttem expansion and there is no 
signature to distinguish between the two. 

The paper is organized as follows. In section |^ we summarize the analysis of quarkonium decay 
into photons and hadrons in the framework of NRQCD. Section |^ describes the NLO calculation 
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and the technique used to isolate and cancel/subtract IR and colhnear divergences. In section ^ 
we give estimates for the non-perturbative matrix elements by comparing the NLO predictions for 
total decay rates with experimental data. Finally, we present a numerical study of the impact of 
octet states on the shape of the photon spectrum. The last section is devoted to our conclusions. 
Appendix |^ collects symbols and notation, appendix ^ collects the results for the Born decay rates 
in D dimensions. A summary of the NLO results is provided in appendix 0, where differential 
decay rates are presented in their final form, after cancellation of all singularities. 



2 NRQCD and fragmentation 

A consistent description of the photon energy spectrum in T ^ 7 + X decay requires the inclusion 
of the fragmentation components The differential photon decay can be expressed in terms of 
a convolution between partonic kernels Ca and the fragmentation functions Da^Y- 



dV „ _ ^ /■! dx 
dJ 



— Jz X X 

a=q,q,9 

= C^ + YCa®Da^^, (1) 

a 

where z = E^/rriQ is the rescaled energy of the photon {rriQ is the heavy-quark mass). The first 
term corresponds to what is usually called the 'prompt' or 'direct' photon production where the 
photon is produced directly in the hard interaction while the second one corresponds to the long- 
distance fragmentation process where one of the partons fragments and transfers a fraction of its 
momentum to the photon. 

Each type of parton, a, contributes according to the process-independent parton-to-photon 
fragmentation functions D^_^^ and the sum runs over all partons. Note that although the frag- 
mentation functions are non-perturbative, we can assign a power of coupling constants, based on 
naively counting the couplings necessary to radiate a photon: since the photon couples directly 
to the quark, -Dg^^ is of 0{acm), while we might expect that Dg^.y is of O^aemCts)- An explicit 
calculation at leading order in gives: 



zDq^^{z) = el^zVg^^{z)\ogj^ , (2) 
zD,^^{z) = 0, (3) 



where the log(Q^/A^) in eq. (0) comes from the integration over the transverse momentum of the 
emitted photon and A is a collinear cut-off that reveals the breaking of the perturbative approach 
and can be chosen of the order of Aqcd- The photon fragmentation functions evolve with just 
as the usual hadronic fragmentation functions do, as a result of gluon bremsstrahlung and qq pair 
production. Such evolution can be derived from a set of coupled equations, which are the usual 
Altarelli-Parisi equations but with an added term that takes into account the leading behaviour 
in eq. (0). The main result of the evolution is that Dg^^ acquires a non- vanishing contribution so 
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that all the Da^^y show the typical logarithmic growth of eq. (0) . This leads to using the following 
leading-log approximation (LLA) for the fragmentation functions 0: 

D.^,{z,Q) = ^^Uz), (4) 

where fa{z) are to be extracted from the data. This shows explicitly that in general the deter- 
mination of the spectrum at 0(aem«s'^) requires the knowledge of partonic kernels Ca in eq. (|l|) 
at 0{aj^^^). This observation was first made, in quarkonia decays, by Catani and Hautmann |Q 
who evaluated the effects of fragmentation contributions to the photon energy spectrum within 
the CSM. They found a strong enhancement in the region of small where soft radiation becomes 
dominant. 

In the NRQCD perspective, a heavy-quarkonium state is represented by a superposition of 
infinite QQ pair configurations organized in powers of f ; v = {v )^/^ is the average velocity of the 
heavy quark in the quarkonium rest frame. Within this framework, the decay width is expanded 
in terms of the matrix elements of 4-fermion operators (that create and annihilate a given QQ 
pair) times perturbative coefficients associated to each operator. By implementing the NRQCD 
factorization formalism within the fragmentation picture, the effects of higher Fock components 
in the quarkonium state can therefore be evaluated systematically. 

The NRQCD expansion for the coefficients Cj(x) reads: 

C^ = T.C^[Q] 1 = 1A.1.9. (5) 

Q 

cm = a[Q](«s(mQ),/iA)i^^H§^^, (6) 



m 



where fi\ is the NRQCD factorization scale and C'j[Q](a;, asimq), the perturbative coefficients 
(here we have dropped the dependence of C*j on the fragmentation scale /xp)- The NRQCD sum 
is performed over all the relevant spin, angular momentum and colour configurations Q that 
contribute at a given order in v. In the case of a T, the structure of the Fock state at order v'^ is 

|T) = 0(l)|66p5fl]) + J2 Oiv)\bb{'PP]) + 0{v^)\bb[^S^^^]) + 0(i;2)|66p5i''"] ) • (7) 

J 

As a consequence, eq. (^ can be written in the following explicit form: 

* * * J * 

+ C^l^S^^^] ^'^^^'^^^'^^^^ (8) 

Let us consider the direct contributions (2 = 7). The leading colour-singlet dimension-6 operator 
contribution is of O(a^aem), and the Vi — operator contribution is suppressed by v'^. All the 
colour-octet processes start contributing at 0{asaemV'^)- By naive power counting, and using the 
approximate relation ~ v"^, one finds therefore that the octet states contribute to the same 
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Figure 1: Sample of LO Feynman diagrams: direct and fragmentation. 



order as the singlet relativistic corrections and might be comparable in size to these. Moreover 
differential quantities are obviously sensitive to the details of the kinematics and so it may happen 
that contributions that are suppressed by standard counting rules are actually leading, in some 
particular region of the phase space. 

LO diagrams are shown in fig. |l|. By considering the following perturbative QCD expansions 
of the coefficients Ca[Q] and of the fragmentation functions Di^j, 



A^, = A^°!, + A^H, + ---, (10) 

one is able to write the general structure of the LO spectrum: 

dr(o) 



dz 



{Cf [Q] + [Q] ® + 2Cf [Q] ® Djo)^} . (11) 



Since the LO colour-octet contributions have a two particle final state, the kinematics is fixed and 
the delta function 6{1 — x) of the short-distance coefficient transforms the convolutions in trivial 
products: 



d^ 



E [rBorn(Q gi)S{l -z)+2 FBorn(Q ^ 99)D^'i-y{z) 

Q 

+2j2Tsorn('S? ^qq)D^^i^{z), (12) 



where the first sum is performed over the lowest-order non-zero octet configurations Q = ^5*0 , 
^Po'^'; ^-^2^^^ while the second one over the flavours of the light quarks. As eq. (O) shows, at 



leading order the colour-octet contributions are proportional to the fragmentation functions and 
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Figure 2: Fragmentation functions of a light parton into a photon according to the 
reference j^. 



to terms proportional to S{1 — z) which do not contribute for z < iQ. 

The fragmentation functions of a hght parton into a photon have been calculated by several 
groups P, H]. In this paper we employ the set recently developed by Bourhis, Fontannaz and 
Guillet 0. In fig. ^ functions Dg^^{z) and J2qDg^^{z) are shown: as was previously stated, the 
contribution from quarks is dominant. 

3 NLO radiative decays: the calculation technique 

In this section we briefly describe the strategy for the calculation of higher-order corrections. A 
consistent calculation of these entails the evaluation of the real and virtual emission diagrams, 
carried out in D dimensions. The UV divergences present in the virtual diagrams are removed by 
the standard renormalization. The IR divergences appearing after the integration over the phase 

^Although we did not include these 'direct' terms in our analysis, we expect that resummation of higher order 
effects for z ^ 1 will induce an effective smearing of the delta function and 'feed down' some photons to lower 
values of z p|. This point will be discussed in more detail in the sequel. 
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space of the emitted parton are cancelled by similar divergences present in the virtual corrections, 
or by higher-order corrections to the long-distance matrix elements [Q. Collinear divergences, 
finally, are either cancelled by similar divergences in the virtual corrections or by factorization 
into the NLO fragmentation functions. The evaluation of the real emission matrix elements in D 
dimensions being particularly complex, we follow in this paper the technique developed in ref. [|T^ 



and already employed in |jTT|, |I2|, whereby the structure of soft and collinear singularities in D 
dimensions is extracted by using universal factorization properties of the amplitudes. Thanks 
to these factorization properties, the residues of all IR and collinear poles in D dimensions can 
be obtained without an explicit calculation of the full D-dimensional real matrix elements. In 
general they only require the knowledge of the Z)-dimensional Born-level amplitudes, a much 
simpler task. The isolation of these residues allows the complete cancellations of the relative 
poles in D dimensions to be carried out, leaving residual finite expressions, which can then be 
evaluated exactly directly in D = 4 dimensions. In this way one can avoid the calculation of the 
full D-dimensional real-emission matrix elements. Furthermore, the four- dimensional real matrix 



elements that will be required have been known in the literature for quite some time ^3], |T^. The 
study of the soft behaviour of the real-emission amplitudes was already presented in |[Tl|, |12| and 
we made substantial use of those results. 

To be more specific, let us consider the three-body decay processes ki + k2 + k^, where 

= QQ[^^~^^Lj^^'^\ Using the conservation of energy-momentum and rotational invariance, 
it is straightforward to verify that there are only two independent variables, which we chose to be 
Xi, the fraction of energy of the parton whose spectrum we are interested in, and y , the cosine 
of the angle of such parton with one of the other two. Within this choice, the differential decay 
width in D dimensions reads : 

' -^MQ] 6(1 - X,) ^ )[Q] + Cf . (13) 



2MSo 



The NLO spectrum coefficients are the sum of the virtual and the real [R) and the virtual {V) 
QCD corrections. In general both channels ggg and qqg contribute to the real term, the Si^2 are 
factors that account for the right counting for identical particles in the final state, and for the 
multiplicity of the various corrections, and $(2) is the total two-body phase space in D dimensions: 



1 f^T^y T(l-e 

while N and K are defined as 



^(2) = -(t7^) (14) 



n\M^J r(2-2e 



M ^ / 47r \ TT^ 
iV=(^(^) r(l + e), i^ = r(l + e)r(l-e)~l + e^-^. (15) 

The function f{x,y) is defined as 

fR[Q]{x^,y) = il-xM^~y)T.\MQ]i^^,y)f (i6) 

fv[Q] = 2Re£ (An A*y). (17) 
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Figure 3: Sample of NLO Feynman diagrams: direct and fragmentation. 



7 



Since divergences can appear only at the border of phase space, i.e. y = 0,y = l,Xi = 0,Xi = 1, 
ffi is finite for all values of x and y within the integration domain. Therefore all singularities of 
the total decay rates can be easily extracted by isolating the e — poles from those factors in 
eq. (|13D that explicitly depend on Xi and y. It must be noted that an infrared divergence arises 
in the limit Xi —>■ when i = g, giving a term of the form ~ loga;^ in the width. Nevertheless we 
are not interested in regularizing such a divergence, since, in this case, the physical resolution of 
the detector works as a physical cut-off. For the same reason the virtual gluon emission at Xg = 
has not been included in the account of the multiplicities. 

The virtual coefficients can be extracted straightforwardly from ref. |T^. The calculation of 
the real coefficients is much more complicated, and it has been carried out by exploiting the soft 
properties of the amplitude obtained in refs. |1TT| , |12| . To illustrate the fundamental steps of the 
calculation of the real part, we consider here the C^^''[Q] coefficients, with Q being one among 

the C-even configurations ^Sq'\ ^Pq'\ '^-P^^^- Let also Uf = for the time being, so that we neglect 
contributions coming from the decay into qqg. In this case we reorganize the first term of eq. (p!7| ) 
by expanding the structure in powers of e and using the symmetry of the phase space. Considering 
the spectrum of the gluon "1" , we find 



${2) 1 



+ 2x 



2M K Si 

log(l — X 



1 



/i?[Q](x,0) (-— + 2\ogx) 



1 — X 



1 — X / + V ecou 

fR[Q] (x, 0) - -5(1 - x) [' dy[y{l - y)]-'-'fR[Q] (x, y) 
e Jo 



1-xJ + Jo \y 



The soft divergences ~ 5(1 — x) cancel by adding the virtual contribution in the same area 
of the phase space. The last piece of eq. (|1^) is a state-dependent finite contribution. The limit 
y corresponds to gluon 1 and gluon 2 becoming collinear 1||2 and the factor 2 in front 
accounts for the case 1||3. Integration over the phase space gives rise to a pole labelled by Ccou 
and a universal finite part. This divergence is not cancelled by adding the virtual term and 
reveals that non-perturbative effects are leading in this case. In fact the residual sensitivity can 
be consistently factorized into the fragmentation function of the gluon into the photon. Such 
singular residual collinear part corresponds to the first term in eq. (^) plus the collinear piece of 
the virtual contribution (~ 6{1 — x)) that comes from the gluon, ghost self-energy loops of the 
gluon we are selecting, so that it reads 



1 



Ccoll 

2Ca 



r(l + e)^x 

TT 



X 



X 



'1-x) 



+ x(l — x) 



X 
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Ca5(1-x) rBor„[Q]. (19) 



If we now switch on the light fiavours including qqg and gluon vacuum polarization diagrams, then 
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Figure 4: Different colour-octet contributions to the photon spectrum in the T decay up 
to 0{v'^). The differential decay widths dF/dz are reported as a function of z = E^/m. 
All the distributions displayed are normalized to the respective Born (that is Q — s- 57 for 
C-even and ^sf'^ — > qq for the only C-odd contribution). 



we obtain the conventional counter-term ~ Vgg, which has to be subtracted at the factorization 
scale /ip. 

This procedure can be extended to all short- distance terms and may be useful to express the 
factorization in a more general way. At NLO the individual terms in eq. (|I|) may be divergent and 
will be denoted by tilded quantities. As we have already mentioned, such divergences correspond 
only to two final partons becoming coUinear, and their form is dictated by the factorization 
theorem. According to this we can reorganize them as follows: 

C7 = + Cq ® Qq—t-y + Cg ® Qg^-y ; 

g 

Cq = J2 Cq' ® Qq'^q + Cg® Qg^q, 
q' 

Cg = T^Cq^Gq^g + Cg^Gg^g, (20) 

q 

where all of the divergences are now concentrated in the factorization-scale-dependent, transition 
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functions Q 



Ga^b 




Vbaix) 



+ Kab(x) 



(21) 
(22) 
(23) 



where a = g,q,q and all the coefficients Ci are now finite for e —>■ 0. The functions Vba{x) are the 
D = 4 Altarelli-Parisi splitting kernels, collected in appendix 0, and the factors Kij are arbitrary 
functions, defining the factorization scheme. In this paper we adopt the MS factorization, in which 
Kij{x) = for all The coUinear factors Qi^j are usually absorbed into the bare fragmentation 
functions by defining 



D 



'9^7 ^ / y yq^q' 
q' 



d3 



B 



D 



D 



B 

9^1 ' 



(24) 



so that we can write the physical decay rate in terms of finite quantities, 

dr 

q 



(7 + X) = + ^ ® Dq^^ + Cg®D 



As illustrated in fig. 

dr(i) 



1, we write the general structure of the NLO processes as 



dz 



Q 

+ 2C«[Q] ® + 2Cf 



(1) 
9^7 



(25) 



(26) 



4 Results 



In the previous section we have shown how the short-distance coefficients have been calculated 
and all the final-state coUinear divergences have been consistently absorbed into fragmentation 
functions. Now, in order to investigate the phenomenological applications of colour-octet states, 
an estimate of the NRQCD matrix elements (ME) must be given. The long-distance MEs can be 
calculated on the lattice, extracted from experiments when enough data are available, or roughly 
determined by using scaling rules of NRQCD or by renormalization group (RG) arguments. At 
the present time none of the aforementioned techniques is able to give a set of precise values for 
MEs and, as we will see below, estimates are affected by large uncertainties. 

The velocity-scaling of the MEs is basically determined by the number of derivatives in the 
respective operators and by the number of electric or magnetic dipole transitions between the 
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Figure 5: Ratio r(T had)/r(T — > /u"^^") versus renormalization scale fi for differ- 
ent values of A5. The solid lines include NLO colour-octet contributions with the RG 
estimates of the matrix elements. The dotted lines include colour-singlet only. 



QQ pair annihilated at short distance and the QQ pair in the asymptotic physical state. This 
can nicely be described by a multipole expansion of the non-perturbative transition T — >■ Q: an 
^5*0^' can be reached by a chromo- magnetic dipole transition, an ^5*1^^ by a double chromo-electric 
emission and ^-Pf' by a simple chromo-electric transition. The first two are of order while the 
last only of order v^. Finally, since the hard-production vertex for a P-wave is already suppressed 
by f ^ relative to the production of an 5'-state, one realizes that the colour-octet C-even states and 
^S^^ all contribute at the same order in v. Following this approach, we can write: 

(T|C>8(^^i)|T) ^ v"" {T\OiCSi)\T) {T\OsCPj)\T) ^ mV {T\OiCS,)\T) , (27) 

where for bottomonium one usually takes t;^ ~ 0.1 and m ~ 4.8 GeV. 



An alternative approach has been considered by Gremm and Kapustin in |]T5|. They obtain 
estimates for the colour-octet operators by solving the RG equations. To order f ^ and leading 
order in a^, they read: 

A^{T\OsCSo)\r) = Oia^v'), (28) 
aA 
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Figure 6: Ratio r(T — > had)/r(T — > fi^fj-^) versus colour-octet matrix elements for 
different values of A5. The dashed lines indicate the 2a interval of the experimental value 
for Rfj,. 



A^{T\OseS,)\T) = ^l^(T|08fPo)|T), (29) 

(aj1\. TXTTXi 



A^{r\OsCPo)\r) = ^^imvr{T\0,i'S,)\T), (30) 
aA Sivr 

where we used the heavy quark spin symmetry to reexpress the expectation values of 08(^-^1,2) in 
terms of Os{^Po)- We note here that our normaUzation for the colour-singlet NRQCD operators 
differs from the original one introduced by BBL, i.e. Oi = ^^r^f^^- Equation (^) differs from 



the respective equations that appear in ref. [|T5] because we included the contribution of ^P|^' to the 



evolution of Os{^Si), which was left out in the previous treatment 0. Assuming that logarithmic 
terms of the evolution are dominant over the MEs evaluated at a starting scale A ~ Aqcd, 
we obtain: 



^ The authors of ref. agree that it is correct to include the ^p|^' contribution in the right-hand side of eq. 



( |29| ) (private communication). 
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Figure 7: Various Fock contributions to the photon spectrum as a function of z = E^/m. 
The sohd hne gives the LO singlet contribution. Fragmentation and NLO direct are 
summed up for each colour-octet state. The NRQCD MEs are related to the colour-singlet 
one through the RG estimate. The colour-singlet matrix element is arbitrarily chosen to 
be (T|Oi(^S'i)|T) = M^/47r, so that comparison with ref. ||5[ is straightforward. 



« ^"'£'°<sJ;;^))^^|o.('50IT). (31) 

(T|Os(»Po)|T),, . ^,„=„.^,og(-i-^)(T|O,(3S0|T>. (32) 

{T|08('S„)|T)^o « O- (33) 

Once numbers are plugged into the previous expression, one realizes that MEs in eqs. (|27|) 
result larger by more than one order of magnitude with respect to the RG estimates shown in 
eqs. (pID-(p3D. This suggests that the very first assumption, i.e. that the non-perturbative matrix 
elements should be dominated by QCD evolution, is doubtful and cannot be justified unless their 
input values were accidentally much smaller than the 'natural' values given in eqs. (]27|). In any 
case, therefore, the estimates in eqs. (PT|)-(|5BD provide a lower limit for the range of all possible 
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values. 

To obtain an independent test of the RG estimates and possibly find an upper limit for the 
matrix elements, we have analysed their impact on the observable -R^(T) = r(T had)/r(T 
fi~^fi~). This observable is particularly advantageous because of the cancellation of several sources 
of uncertainties: both the colour- singlet NRQCD matrix element and the overall dependence on 
the bottom mass cancel in the ratio. As a result the mass enters only in the logarithm of the 
renormalization scale and its uncertainties can be naturally associated to the choice of the scale 
itself. 

In fig. ^, the ratio is plotted versus the renormalization scale fi (i.e. the NRQCD factoriza- 
tion scale is kept equal to the renormalization one); T^^ includes the 'Pi(^S'i)-operator contribu- 



tion and the NLO QCD corrections |jT6|; Fhad includes the NLO QCD colour-singlet and the 



Pi(^S'i)-operator contribution |T5| (dotted curves) and the NLO colour-octet ^ Sf\^ S^o\'^ P^j 



J 

[|lT] (solid curves). The dashes lines limit the 2a band of the experimental value of i?^ = 37.3 ±1.0 
[|19]. The theoretical curves are drawn according to the following choice of parameters: = 0.1 
and aEui^b) = 1/132. Hence fig. ^ shows that, once the colour-octet RG MEs estimation is 
plugged in, the ratio is consistent with the experiments only for A5 ~ 140 MeV {as{Mz) — 
0.110). On the other hand if we drop the colour-octet term, just the NLO colour-singlet contribu- 
tion can still reproduce the experimental measure of i?^ by choosing a much higher value of A5, 
namely A5 ~ 220 MeV {as{Mz) ~ 0.118). 

Now we fix the renormalization scale /xr = 10 GeV. We note that, more than corresponding 
to the 'natural' choice /ir ~ Mx, this value also satisfies the so-called 'minimal sensitivity princi- 
ple' [^, i.e. it is the value at which /iR^^^i?^(/iR) vanishes. Within this choice, we plot the ratio 



R^ versus the variable 

^ (T|08(3Po)|T) ^ {T\OseS^)\T) 

The result is shown in fig. The solid lines represent the theoretical calculation of i?^ and the 
dashed lines are the 2a experimental range, as in fig. ^. The larger the colour-octet MEs are, the 
smaller A5 has to be taken. In particular, already for values of the MEs of the order of twice the 
RG estimates, we would find a value of A5 ~ 80 MeV {as{Mz) — 0.102), well outside the present 
world average range. 

Following this line one finds that the MEs provided by the velocity scaling rules are strictly 
excluded. In an ideal global fit perspective both the value of A5 and the colour-octet MEs should 
be extracted from the data. Unfortunately the experimental inputs in the T decay sector are not 
sufficient to perform a fit of such a large number of unknown parameters. 

As a confirmation of what we found in fig. |^, fig. ^ shows that the RG estimate reproduces the 
experimental value of i?^ for A5 ~ 140 MeV. Such a value of A5 corresponds to as{mb) — 0.190 
and as{Mz) ~ 0.110. The world average of a, (as(M^) = 0.119±0.004) (or equivalently A5 ^ 237 
MeV) is actually consistent with a vanishing (or even negative) octet contribution to the T decay 
into hadrons. Nevertheless the uncertainties involved are still large: NNLO QCD corrections 
(reflected in the n dependence of the NLO correction) might be important as well as higher twist 
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effects. A clear indication that higher order effects are not neghgible, comes from the two-loop 
calculation of the leptonic width recently performed by Beneke et al. in this case, it is found 
that the 0{a'l) corrections (NNLO) are of the same size (or even larger) of the NLO ones. 

Summing up, we can say that on the one hand comparison with scaling rules of NRQCD shows 
that RG estimates have to be thought of as a lower limit, while on the other hand consistency 
between theory and experiment in total decay rates strongly disfavour much larger colour-octet 
MEs. We then conclude that the RG estimates of the colour-octet MEs, although sizebly smaller 
then expected from NRQCD scaling rules, are the most reasonable at the present stage of our 
knowledge. 

In fig. ^ we show in detail the contribution of the single colour-octet components. The figure 
reports LO, direct, and full NLO contributions for states normalized to their respective Born 
decay widths at 0{as<yem)- Let us consider the C-even states first Csj^\^PS^\^Pl^^). It is evident 
that they contribute to the spectrum with a very similar shape: there is a strong enhancement 
at low values of z due to the fragmentation contribution that is present both at LO and NLO. 
Then it is clearly seen that direct photons mainly contribute near the end-point, a zone of the 
spectrum where the fixed-order calculation is not reliable: in fact there are clear indications of 
a need of resummation both in the short-distance perturbative expansion in a.^ and in the long- 
distance V series. In ref. Rothstein and Wise identified an infinite class of NRQCD operators, 
which determine the shape of photonic end-point functions, and introduced the so-called 'shape 
function', to be extract from data. The overall effect of colour-octet states would be a smearing 
of the energy distribution near the end-point on the interval ~ 0.1. In the case of the 
component, the direct amplitude is not divergent in 2; = 1 and the NLO correction to the LO 
fragmentation picture is very small. Indeed the NLO contribution from direct photons is negative 
in the MS-renormalization scheme and is almost balanced by the other NLO fragmentation terms. 

Finally figs. and |^ show the total contribution to the spectrum, using the RG estimate for the 
non-relativistic matrix elements. We notice that the overall effect of octet states is at its minimum 
in the central region of the spectrum, exactly where the singlet LO direct contribution dominates. 
This indicates that this region of the spectrum is 'safe' from colour-octet effects, and therefore 
we think that it should be used to make a comparison with experimental data. Moreover this 
indicates that relativistic corrections to the singlet (which are indeed important) and higher-order 
strong ones should be included to have a consistent theoretical picture at NLO. On the other side, 
for small values of z, colour-octet components are not negligible. In this area of the phase space, 
the fragmentation components from gluons contribute at the same order in as the ones from 
quarks, and there is no signature to distinguish between the two. Contrary to LO expectations in 
the framework of CSM 0], we conclude that the decay of T into a photon would not be useful for 
an estimate of the photon fragmentation functions. 

As a final remark, we notice that, not surprisingly, many of the aspects of the photon spectrum 
in quarkonium decay, resemble those in photoproduction p2|, |2^. Cross-sections plotted versus 



the inelasticity z of the quarkonium state show a very similar pattern: for z ^ 1, a divergence, 
which is not supported by the available experimental data, reveals the breaking of the NRQCD 
expansion in powers of as and v. On the other side, for low values of z, the resolved contributions. 
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0,2 0.4 0,6 0,8 

z 

Figure 8: Total colour-octet contribution on the LO, colour-singlet photon spectrum. 
Notice that neither NLO QCD nor relativistic effects are included in the singlet contri- 
bution. Normalization and MEs as in fig. 0. 



which corresponds to fragmentation in the decays, are indeed dominated by colour-octet states. 

5 Conclusions 

We presented the calculation of C(a^aem) colour-octet corrections to the decay of T into one 
photon plus light hadrons. Both direct and fragmentation contributions have been included at 
NLO. In order to study the impact of these contributions on the photon spectrum, an estimate 
of the non-perturbative MEs was also given. By comparing the available experimental data on 
fully inclusive and leptonic decay rates with the NLO theoretical predictions of NRQCD, we 
found an unexpected result: estimates based on naive scaling rules result in large colour-octet 
contributions to the total rates which are not consistent with the data. In particular, it turns out 
that non-perturbative MEs should be much smaller then expected from NQRCD scaling rules. 
Nevertheless, using the above mentioned estimates for the non-perturbative MEs, we showed that 
there are sizeable effects at the end-points of the spectrum of the photon. In the case of low values 
of z, the possibility of measuring the fragmentation function of a gluon into a photon, which was 
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suggested by the LO result in the CSM p, becomes unfeasible: for the colour-octet states both 
quark and gluon fragmentation processes are of the same order in agOteui and there is no signature 
to distinguish between the two. Moreover, for values of z near the end-point, breaking of the 
fixed-order calculation is manifest, and the resummations of both short- distance coefficient in 
and non-perturbative MEs in f , are called for. Nevertheless a 'safe' region, for 0.3 < z < 0.9, has 
been found where octet effects are at their minimum and the perturbative expansion in powers of 
Us and V under proper control. Following this point of view, we consider the NLO QCD correction 
the colour- singlet differential decay dr/d-E'^(T — » ^^i'^' 199) worth while to be undertaken. 
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A Symbols and notations 



This appendix collects the meaning of various symbols, which are used throughout the paper 
Kinematical factors : 

M = 2m , V = \ I 



M2 



(35) 



where s is the partonic centre-of-mass energy squared and Shad is the hadronic one; v is the velocity 
of the bound (anti) quark in the quarkonium rest frame, 2v then being the relative velocity of the 
quark and the antiquark. The following expression is used: 



r(l + e) = 1 + e i-^E + log(47r) + log ^ j + 0{e^) . 



(36) 



Altarelli-Parisi splitting functions. Several functions related to the AP splitting kernels enter in 
our calculations. We collect here our definitions: 

"l + a;2 



Pqq (x) 
Vqq{x) 



el-x) 



1 — X 

+ (1 - x)^ - 2e x{l - x) 



Tf 

Tf ^ 
l + (l-x)^ 

X 



x^ + (1 -a;)^ 



— ex. 



(37) 

(38) 

(39) 
(40) 

(41) 
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X 



P,,{x) = Cf 

P,,{x) = 2Ca 
V,,{x) = 2Ca 



1 + (1 -x)^ 



— ex 



X 



1 + (1 -x) 

X 

X 1- X . 
+ + x{l-x) 

1 — X X 

X 1 X . . 

7" ^ H + x{l - x) 

(1 — x)+ X 



+ ho5{l-x). 



(42) 
(43) 

(44) 
(45) 
(46) 



The Pij are the D-dimensional sphtting functions that appear in the factorization of colhnear 
singularities from real emission, while the functions Vij are the four-dimensional AP kernels, 
which enter in the MS colhnear counter-terms. The and 'a' distributions are defined by: 



dx [r(x)]^0(x) = dxr(x) [0(x) -0(1)] , 
/' dx [T{xX<P{x) = [' dx T{x) [<P{x) - 0(1)] , 

J a J a 



(47) 
(48) 



where T{x) is the function associated to the distributions [T'(x)]+,a. We recall a useful weak 
distributional identity: 



[T{x)]^ = [T{x)]^-5{l-x) rT{x)dx. 

J 



In particular it is straightforward to get: 



,1 — a; 
'log(l — x) 
1 — X 



1 X/ a 

logfl — X 



+ 5{l-x) log(l-a), 

1 



+ 



+5(l-x)-log^(l-a) 



(49) 



(50) 
(51) 



Colour coefficients 



2Nr 



Ca = N, 



The following standard symbol is used: 



11 2 



Tf = - 
2 



(52) 



(53) 



with Uf the number of flavours lighter than the bound one. 
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NRQCD operators. To denote a perturbative QQ state with generic spin and angular momentum 
quantum numbers, and in a colour-singlet or colour-octet state, we use the symbol: 



Notice that, according to the discussion in ref. ]TT|, our conventions differ from the Bodwin, 
Braaten and Lepage ones Q (labelled here as BBL) in the case of a colour-singlet. 



2N, 



BBL 



(55) 
(56) 



B Summary of lowest order results 

B.l Born widths 

The decay rates read 

r(T ^ ^ ab) = f(Q[i'«l a6)(T|0[i,8](25+iLj)|T) 



(57) 



the short-distance coefficients T having been calculated according to the rules of ref. [O] . We shall 
use the short-hand notation 



T(Q 



1,8] 



ab) = T{T^Q 



[1,8] 



ab) 



(51 



to indicate the decay of the physical quarkonium state H through the intermediate QQ state 



g[i.8] = QQ[2^^+iLj' '*^]. The D-dimensional (D = A - 2e) 0(a,aem) level decay rates read: 

'^".^^ <I>(2)(1 - e)(l - 2e){T\OsCSo)\T) , 



Born 




8] _ 


9l) 


Born 


es. 


8] _ 


9l) 


rBorn 




SI _ 


9l) 


rBorn 




SI _ 


9l) 


l^Born 


CP2 


SI _ 


9l) 



288asaem e?, u^tt"^ 1 - e , , , , , 

<|.(2);^^(T|08(3Po)|T), 



'3-2e 
0, 

64a,aem /x^^vr^ ^ (6 - 13e + 



(2) 



(3-2e)(5-2e) 



{T\OsCP2)\T). 



(59) 
(60) 
(61) 
(62) 
(63) 



Lowest Ofa?) contributions: 



rBorn('5o''^ ^ ^7^?) = ^^"^l^, $(2)(1 " e)(l - 2e) {T \OsC So)\T) 

(T|08(35i)|T), 



J- Bornl '-'1 



'3-2e 



(64) 
(65) 
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TBorn 




gg) = Bp 


TBorn 




g9) = o, 


TBorn 




gg) = Bp 



lUaln^^7i\ (1-e 

— $^9^-^ 

m4 (3 - 2e) 



*(2)^^^(T|08ePo)|T), (66) 

(67) 



32a^, fj.^'^n^ (6 - 13e + 4e' 
(3-2e)(5-2e 



'^(2) ;, .,,7 .; (T|08(^P2)|T) , (6^ 

where $(2) is defined according to eq. ([I^) . 

B.2 The LO spectrum coefficients C^^^[Q] 

We can now read out the lowest-order coefficients according to eqs. (^-(|ll|). For Q = , ^5*0^^ 
we have: 

Cf[Q](z) = T^ornlQ ^ gim - z) , (69) 

Cf[Q](x) = 2Tj,,,^[Q^gg]6{l-x), (70) 

Cf[Q](a:) = 0, (71) 



and for ^S^^ 



CfpM'^K^) = 0, (72) 
Ci'^?Sf^]ix) = 0, (73) 

Cfp5fi](x) = rBorn[Q-gg]5(i-x). (74) 



C Summary of O(ajaem) results 
C.l The NLO photonic coefficients C^fQ] 

We summarize the NLO spectrum coefficient following the convention of eqs. (|5|)-(p6D.The photon 
energy fraction is 2; = E^/m. Components ~ 6{z) have been neglected. For Q = ^Sq'\ ^P(t\ ^P2^\ 
we have: 



271 



(75) 



where 



Afsf) = (-10 + + (i§ - ^) - Bn,T, . (76) 



20 



14 



Cf -- + 



+ Ca 



'85 TT^' 



10 



18 2 1 9 ^-^"^^ ' 



-8C,. + CA(y +log2 



47 



45 



and 



[8].,. ^ (-72 + 144z-176z^ + 104z3-23^^) , 2 

fiiS'o %z) = Ca ^ , ,2 _ + ^fTp^z , 



6{-2 + zfz 

^ (-960 + 3360 z - 6224 z'^ + 5312 - 1544 z^ - 520 , 

2 



(77) 

(78) 

(79) 



54:{-2 + z)^z^ 
+mz^ - 136/ + 9^^) +r^/^F^(^ + 2)^ 

il^Pf\ {z) = Ca — (-240 + 1848 z - 7820 z'^ + 13976 z^ - 12710 z^ + 6254 

36 (—2 + z) z^ ^ 

-1628 + 197 / - 15 + UfTp— (lO - 5z + z'^) , 

^7ric[8]u 2 (+12 - 36^ + 56;^^ _ 52;^3 ^ 28^^ - 8^^ + z^) 

n^pt\z)^CA 



-2^zY z- 



3 .2 



9 (-2 + zf z^ 
+104 - 66 / + 40 z^ - 10 + , 

7r3p[8]^' . „ 1 



(+160 - 720 ^ + 1624 z^ - 2016 z^ + 1360 / - 468 , 



r2?pr]{z)^CA-, 



3 {-2 + zf z" 
+1826 / - 582 z'' + 115 - 13 + . 

For the component we get: 

1 



(+40 - 348 ^ + 1618 z^ - 3684 z^ + 4702 - 3669 z^ 



(80) 

^5 

(81) 
(82) 

(83) 
(84) 



[8]l 



20Q;emeQQ;f 



+ 



-1 + ^)(8 - 12^ + hz^) log(l - z) 



(T|08P-5i]|T) 



(-2 + ^)3^2' 

+ V^o.A^Sf^qqp^P^,{z) (log^ + log(l-;.) + 2 1og;.We^, 



(85) 



and finally for ^Pf: 

(-7(l)[3p|8]j ^ 2Q;emeQQi3 



1 



- 276z^ + 3l/) + 



-2 + ^)4^3 

12 



{-2^zfz 



(240 + 312z - 2620^^ + 4204^^ - 3150z^ + 1260^^ 
'-\ + ^)(40 + 52^ - 430^2 + 716^^ - 588^^ 



21 



+ 275^^ - 74^6 + Uz'^ - z^) log(l - z) + -Uf- — - 

3 X 



(T|e»8pPi]|T) 



(86) 



C.2 The NLO gluonic coefficients C^^fQ] 

In this section we present the NLO QCD spectrum of the gluon arising from the colour-octet 
components. Contributions ~ 5{x) have been neglected. The gluon energy fraction is denoted by 
X = Eg/m . For Q = ^-Sf ^Pj^', =^Pj^l, we have: 



IT 



4m^^ , , ^, ^ / X /log(l — 
log -r^ssi^) + 2 ^ogxVggix) + 7 ^ 



(1-X)P,,(X) 



+ 



1 — X 



f[Q]ix) + (P[Q] + Abo log ^) 5(1 - x) 



where: 
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pppf] 



and furthermore 



= -8Ci. + Ca (5 + ^ log 2 - iyr^) -^n/T^, 



— (-120 + 336 a; - 494 a;^ + 410 a;^ - 215 a;^ + 72 x^ - 12 



(87) 

(88) 
(89) 
(90) 



X 



+ 



6{-2 + x) X 
2 Ca (—1 + x) 

/ ^ o ^(16 - 40 a; + 50a;2 - 26a;^ - 8a;^ + 16a;^ - 7a;*^ + x'^) log(l - x) 

(2 — x) a;^ 

(91) 



+ nfTp-x, 

Ca 



-4 — (-1536 + 5376 x - 9632 + 10016 x^ - 9288 + 12976 x^ 



+ 



5A{-2 + xyx^ 
16906 + 13918 - 6623 a;^ + 1664 - 172 a;^°) 
2 C A ( — 1 + x) 

^ ^ (256 - 896 a; + 1504 x'^ - 1008 x^ - 516 x^ + 1792 a:^ - 2276 a:^ 

9(2-a;) x4 

2011 a;^ - 1220 x^ + 464 x^ - 99 a;^° + 9 a;^^) log(l - x) 
2{2 + xy 



+ UfTp 



27x 



(92) 
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-384 + 3072 X - 12704 + 25376 x-^ - 30738 x^ + 26998 x^ 



+ 



36(-2 + x) x3 
19231 + 10924 - 4373 + 1028 - 106 

Ca{-1 + x\^^^ -5i2x + 2032x2 - 3840 + 3747 - 1577 - 515 



3(2-a;)^a;4 

+ 1162 - 800 x^ + 308 x^ - 66 + 6x^^) log(l - x) 
+ rifTp-^ (l0-5x + x'^) . 



(93) 



For the component we get: 



4(-l + x) 

2 

[8] 



+ 



(1168 - 3264x + 3740^2 - 2200x3 + 693x^ - 108x^) 

(T|e»8P-5i]|T) 



{-2 + xfx 



(-2 + x)2x 

(584 - 1632x + 1904x^ - 1134^^ + 324x^ - 27x^) log(l - x 



TT 



-3 + 3x — X 



X 



Am 



+ l^log + log(l - X) + 2 logxj Pgq{x) 



(94) 



Finally for ^Pf^: 



5a| 
18 



-2 + x)%3 



(384 + 384x - 4192x2 + 7552x3 - 6446x^ + 2876x^ 



12 



485x^ - 141x^ + 82x^ - lOx^) + r^(-l + ^)(64 + 64x - 688x2 ^ ^280x3 



1181x^ + 626x^ - 195x^ + 36x^ - 4x®) log(l - x) + \nf- — - 



(T|08pA]|T) 



(95) 



C.3 The NLO quark coefficients C^fQ] 



We report in this section the quark energy spectrum m. Q ^ qqg decays. The adimensional energy 
of the quark Eq/m is denoted by x. 



1 c[8ll 



TT 



rBom[^'S'f gg] 



Pgg(x)log— 2- 



2^1 „M , ^ric-lsli 



where 



+ 2 x(l - x)7V + Vq,{x) log[x^(l - x)] + j^St\x) 



/,P5n(x)=x(l-x)(l + log(l-x)). 



(96) 



(97) 
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We have 



[8] 



TT 



1^ , 

-Vgq{x) log— 2- 



+ Ci 



+ 



1 + 



1 



1 -X 



f,['S?]{x) + A[^Sr]S{l-x 



[8]l 



where 



^PS! '1 = - + - j + C. + J log 2 - - j - -„,T, + 2 log 1^ 



and 



= C^F J (-4 + x)- Ca^ (5 - 5x + 2x2) + Ca{1 - x){-2 + x) log(l - 



a; 



[3 p[8] 



a; 



(T|af"P;]|T) 



TT 



+ 7^,3 (x) log 

PF 



2a;(l - a;)T^ + \og[x\l - x)]V,g{x) + 



We also have 



where 



and finally 



CW[3pfl]=i?^«35(l- 



log h Oi 

9 ^ 2m 



(T|e»8PPi]|T) 



1 — X 



(T|08pPi]|T) 



27 
2 

9 



ao-g, ai--, a^--, 



x(33 - 72 X + 43 x^) - 3(1 - x) (4 - 9x + 9x^) log(l - x) 



x{3 + 6x- 5^2) + 3(1 - x) log(l - x) 



36 



x(57 - 9Qx + 53x2) _ 3^;^ _ ^)(5 _ i2x + 12x^) log(l - x) 
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